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Abstract. We apply the method of nonlinear steepest descent to compute the 
long-time asymptotics of solutions of the Korteweg— de Vries equation which 
are decaying perturbations of a quasi-periodic finite-gap background solution. 
We compute a nonlinear dispersion relation and show that the x/t plane splits 
into g + 1 soliton regions which are interlaced by g + 1 oscillatory regions, 
where g + 1 is the number of spectral gaps. 

In the soliton regions the solution is asymptotically given by a number of 
solitons travelling on top of finite-gap solutions which are in the same isospec- 
tral class as the background solution. In the oscillatory region the solution 
can be described by a modulated finite-gap solution plus a decaying dispersive 
tail. The modulation is given by a phase transition on the isospectral torus 
and is, together with the dispersive tail, explicitly characterized in terms of 
Abelian integrals on the underlying hyperelliptic curve. 



1. Introduction 

Consider the Korteweg-de Vries (KdV) equation 

(1.1) V t (x,t)=6V(x,t)V x {x,t)-V xxx (x,t), (i,i)elxl, 

where the subscripts denote differentiation with respect to the corresponding vari- 
ables. 

Following the seminal work of Gardner, Green, Kruskal, and Miura [17], one 
can use the inverse scattering transform to establish existence and uniqueness of 
(real-valued) classical solutions for the corresponding initial value problem with 
rapidly decaying initial conditions. We refer to, for instance, the monograph by 
Marchenko [28] ■ Our concern here are the long-time asymptotics of such solutions. 
The classical result is that an arbitrary short-range solution of the above type will 
eventually split into a number of solitons travelling to the right plus a decaying ra- 
diation part travelling to the left. The first numerical evidence for such a behaviour 
was found by Zabusky and Kruskal [37] . The first mathematical results were given 
by Ablowitz and Newell [T|, Manakov [27 , and Sabat [31] ■ First rigorous results for 
the KdV equation were proved by Sabat [31] and Tanaka [33] . Precise asymptotics 
for the radiation part were first formally derived by Zakharov and Manakov [36], 
by Ablowitz and Segur [2], [32], by Buslaev [6] (see also [5]), and later on rigorously 
justified and extended to all orders by Buslaev and Sukhanov [7J. A detailed rig- 
orous proof (not requiring any a priori information on the asymptotic form of the 
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solution) was given by Deift and Zhou [8] based on earlier work of Manakov [27] and 
Its [18] and is now known as the nonlinear steepest descent method for oscillatory 
Riemann-Hilbert problems. For an expository introduction to this method applied 
to the KdV equation we refer to [20] . For further information on the history of this 
problem we refer to the survey by Deift, Its, and Zhou [9]. 

In this paper we want to look at the case of solutions which are asymptotically 
close to a quasi-periodic algebro-geometric finite-gap solution of the KdV equation. 
In this case the underlying inverse scattering transform was developed only recently 
by Grunert, Egorova, and Teschl [T2], [T3], [T3J. So while the initial value problem 
for this class of solutions is well understood, nothing was known about their long- 
time asymptotics even though the first attempts by Kuznetsov and Mikhailov |23] 
date back over 35 years ago. It is the aim of the present paper to fill this gap. In 
case of the discrete analog, the Toda lattice (see e.g. [34]), Kamvissis and Teschl 
[2"T] , [32] (with further extensions by Kriiger and Teschl [55] ) have recently extended 
the nonlinear steepest descent method for Riemann-Hilbert problem deformations 
to Riemann surfaces and used this extension to prove the following result for the 
Toda lattice: 

Let g be the genus of the hyperelliptic curve associated with the unperturbed 
solution. Then, apart from the phenomenon of the solitons travelling on the quasi- 
periodic background, the (n, i)-plane contains g + 2 areas where the perturbed 
solution is close to a finite-gap solution from the same isospectral torus. In between 
there are g + 1 regions where the perturbed solution is asymptotically close to a 
modulated lattice which undergoes a continuous phase transition (in the Jacobian 
variety) and which interpolates between these isospectral solutions. In the special 
case of the free lattice (g = 0) the isospectral torus consists of just one point and 
the known results are recovered. Both the solutions in the isospectral torus and 
the phase transition were explicitly characterized in terms of Abelian integrals on 
the underlying hyperelliptic curve. 

In the present paper we will use this extension for Riemann-Hilbert problems on 
Riemann surfaces to prove an analog result for the KdV equation to be formulated 
in the next section. 



To set the stage we will choose a quasi-periodic algebro-geometric finite-gap 
background solution V q {x, t) of the KdV equation (cf. the next section) plus another 
solution V(x, t) of the KdV equation such that 



for all £ G K. We remark that such solutions exist which can be shown by solving 
the associated Cauchy problem via the inverse scattering transform [12] . 

To fix our background solution V q , let us consider a hyperelliptic Riemann surface 
JC g of genus g £ No with real moduli Eq,E±, . , . , E^ g . Then we choose a Dirichlet 
divisor T>^ x t ^ and introduce 



2. Main results 



(2.1) 




z(p, x, t) 



(2.2) 
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where A E (a Eo ) is Abel's map (for divisors), and S Bo , V_ , and U_ 2 are some 
constants defined in detail in Section [3] below. Then our background solution is 
given in terms of Riemann theta functions (cf. (|3 . 1 3p ) by 



g 

V q (x, t)=E + J3(J5y-i + E 2j - 2 f i j (x, *)) 

(2-3) j f 

= E + J2{E 2j -i + E 2j - 2Aj) - 28 2 x \n9(z( Poo ,x, t)) , 

3 = 1 

where Xj G (E 2 j_ 1 ,E 2 j), j = 1, . . . ,g. 

In order to state our main result, we begin by recalling that the perturbed KdV 
solution V(x,t), x € K, for fixed t G M, is uniquely determined by its scattering 
data, that is, by the right reflection coefficient i? + (A, i), A € <j{H q ), and the eigen- 
values /9fc € R\er(i/q), fc = 1, . . . ,N, together with the corresponding right norming 
constants 7+ > 0, fc = 1, . . . , N. Here 

9-1 

(2-4) <r(ff,) - |J [E 2j ,E 2j+1 ] U [£fc fl> oo) 

denotes the finite-band spectrum of the underlying background Lax operator 

(2.5) H g (t) = -dl + V q (x,t). 

The relation between the energy A of the underlying Lax operator H q and the 
propagation speed at which the corresponding parts of the solutions of the KdV 
equation travel is given by 

(2.6) V ^ = T 
where 

r(A+ie,+) n. 



-I2Rc(ijy 

(2.7) u(A) = lim 



-r, (■ f (A+ie,+) \ ' 

and can be regarded as a nonlinear analog of the classical dispersion relation. Here 
u; Pooj o and w Poci 2 are Abelian differentials of the second kind on the underlying 
Riemann surface defined in (|3.15[) and (|3.16|) . We will show in Section [5] that v is 
a decreasing homeomorphism of K and we will denote its inverse by C( v )- 
Furthermore, we define the limiting KdV solution Vi tV (x, t) via the relation 

/ (V l;V -V q ){y,t)dy = - V 4i / u Poo , + - / log(l - \R\ 2 )u Pee , 
Jx _T7,s JE(pA n JcM 



(2.8) 




E{ Pj ) JC(v) 

(z(poo,x,t) +5(v)) 



with 

S e (v) = -2 MpMPi) + / lo g(l - l^| 2 )C£, 

where i? = i?+(A, t) is the associated reflection coefficient and is a canonical basis 
of holomorphic differentials. Moreover, C(v) is a contour on the Riemann surface 
obtained by taking the part of the spectrum <r{H q ) which is to the left of ((v) and 
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lifting it to the Riemann surface (oriented such that the upper sheet lies to its left). 
Here we have also identified pj with its lift to the upper sheet and E(pj) denotes 
the branch point closest to pj. If v = x/t we set Vi(x, t) = Vi tX u(x,t). 

Then our main result concerning the long-time asymptotics in the soliton region 
is given by the following theorem: 

Theorem 2.1. Assume V{x,t) is a classical solution of the KdV equation (jl.ip 
satisfying 

r+oo 

(2.9) / (l + \x\ 1+n )(\V(x,t)-V q (x,t)\)dx<oo, 



for some integer n > 1 and abbreviate by Ck = v(pk) the velocity of the k 'th soliton. 
Then the asymptotics in the soliton region, {{x, t)\ C,(x/t) G M\a(H q )}, are the 
following: 

Let e > be sufficiently small such that the intervals [c k — e, Ck + e], 1 < k < N, 
are disjoint and lie inside v(M\a(H q )). 

If \j — Ck\ < e for some k, the solution is asymptotically given by a one-soliton 
solution on top of the limiting solution: 

(2.10) J (V-V l , Ck )(y,t)dy = -2—log(c l , k (x > t)) +0(t~ n ), 
as well as 

(2.11) (V-V hCk )(x,t) = 2^\og( Cl , k (x,t))+0(t- n ), 
where 

(2-12) ci tk (x,t) = 1 +7fc / ipi, Ck (p k ,y,t) 2 dy 



an 



d 



Ik = Ik 



(2.13) 



0U(pkA0) + S(ck))Y TT exp (2 ojp . p * 
■expf^i/ log(l - |i?| 2 M 

V m JC(CL.) 



Here ipi^ v (p,x,t) denotes the Baker- Akhiezer function corresponding to the limiting 
KdV solution Vi tV (x,t) and cu pq denotes the Abelian differential of the third kind 
with poles at p and q. 

If \ j — Ck\ > E, for all k, the solution is asymptotically close to this limiting 
solution: 

/■OO 

(2.14) / (V-V l )(y > t)dy = 0(t- n ), 

J X 

as well as 

(2.15) V{x,t) = Vi(x,t) + 0(t- n ). 

In particular, we see that the solution splits into a sum of independent solitons 
where the presence of the other solitons and the radiation part corresponding to 
the continuous spectrum manifests itself in phase shifts given by (|2.13[) . Moreover, 
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observe that in the periodic case considered here one can have a stationary soliton 
(see the discussion in Section [5]) . 

The proof will be given at the end of Section \5\ 

Theorem 2.2. Assume V{x,t) is a classical solution of the KdV equation 
satisfying (|2.ip and let Dj be the sector Dj = {(x, t) : C(x/i) € [E2j + £, Hhj+x ~ s]} 
for some e > 0. Then the asymptotic is given by 



(2.16) £°°(V -V0(y, t)dy = J ^L_R e (/3(x, t))A^zj) + 0(t~ a ), 



respectively 
(2.17) 



V w' 1 1 k=i 1=1 J 

for any 1/2 < a < 1 uniformly in Dj as t — > oo. Here 
(2-18) <t>(j>) = -24i u Poo ,i - 2ij [ cj Poofi 

J Pa 1 Jpo 

is the phase function, 

(219) fM= 12 nu^(^--^) . 

(where -R^g+iM ^ e square root of the underlying Riemann surface IC g and we 
identify Zj with its lift to the upper sheet), 

g g 

(2.20) Af(zj) = u) Poofi (zj) - ^^c w (^)a g _i(^)Cfe(^), 

k=l £=1 

with w Pooi o an Abelian differential of the second kind with a second order pole at 
Poo (cf eq- (|3.15|) ). and uj(j>) denoting the value of a differential evaluated at p in 
the chart given by the canonical projection, and Cki(E), a g -i{0i) some constants 
defined in (|6.24j) . (16.34[) . respectively. Moreover, 

/3{x,t) = y?;e i (' r / 4 - ar s( i? '(^»+ ar s( r ( il/ ))+ 2ra (^)) ^ 7 fr) j ~ 1V c -t4>Ui) t -iv 

6(z(Zj, 0,0)) d(z{z],x,t) + 5[x/t)) 
0(z(zj,x,t)+6(x/t)) 6{z{z* p 0,0)) 

/ fPk If / 1 — |i?| 2 \ 

(2.21) exp - V / lo ZjZ . + — log J ' )uj ZjZ , 

where T(z) is the gamma function, u Zj z * an Abelian differential of the third kind 
defined in (]3.21[) . 

(2.22) v = -±log(l-\R( Zj )\ 2 ) >0, 



and ot(zj) is a constant defined in (|6.9I) . 

The proof of this theorem will be given in Section |6] of this paper. 
Finally, note that if q(x,t) solves the KdV equation, then so does q(— x, — t). 
Therefore it suffices to investigate the case t — > oo. 
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3. Algebro-geometric quasi-periodic finite-gap solutions 

This section presents some well-known facts on the class of algebro-geometric 
quasi-periodic finite-gap solutions, that is the class of stationary solutions of the 
KdV hierarchy, since we want to choose our background solution V q from that class. 
We will use the same notation as in [16 , where we also refer to for proofs. As a 
reference for Riemann surfaces in this context we recommend |15j . 

To set the stage let JC g be the Riemann surface associated with the following 
function 

29 

(3.1) R 1 2 / g li(z)=il[V^E-, E <E 1 <---<E 2g , 

j=o 

where g G No and {Ej}^ Q is a set of real numbers. Here ^J~. denotes the standard 

1 /2 

root with branch cut along (0, oo). We extend R 2g+ i ( z ) to the branch cuts by 

setting i?2g+i( z ) = nm eio R 2 g+i( z + i £ ) f° r z Hence we have 

(3.2) 

r forze (-oo,E ), 
1/2 M _ | ff i/2 , x, I (-1) 9+J i for z€ (E 2j ,E 2j+1 ), j = 0,...,g-l, 
25+11 ' ~ in ^+i( z )\ ■ S for z G (E 2j+1 ,E 2j+2 ), j g 1. 

[ i for z G (E 2g , oo). 

/C ff is a compact, hyperelliptic Riemann surface of genus g. 

1 /2 

A point on K. g is denoted by p — (z,±R 2g+1 (zj) — (z,±), z G C, or = 
(00,00), and the projection onto C U {00} by 7r(p) = z. The points {(Ej, 0),0 < 
j < 2g} U {(00, 00)} C K a are called branch points and the sets 

9-1 

(3.3) n± = {(z,±Ry g 2 +1 {z)) I z G C\ [}[E 2j ,E 2j+l ] U [£29,^)} C K a 

are called upper, lower sheet, respectively. 

Next we will introduce the representatives Wj,bj}^ =1 of a canonical homology 
basis for IC g . For o,j we start near E 2 j-i on n + , surround E 2 j thereby changing to 
n_ and return to our starting point encircling E 2 j-\ again changing sheets. For 
bj we choose a cycle surrounding Eo,E 2 j—i counterclockwise (once) on n + . The 
cycles are chosen such that their intersection matrix reads 

(3.4) dj o cij = bi o bj =0, a, o bj = Si.j, 1 < i,j < g. 

The corresponding canonical basis {Cj'}f=i f° r the space of holomorphic differ- 
entials can be constructed by 

(3-5) O^E^fe)-^, 

fe=l Rlg+\ 

where the constants Cj(k), j, k = 1, . . . ,g are given by 

(3.6) c j{ k) = C ik \ C jk = J ^ = 2/ ^T^M. 

Ja k H 2g+1 JE 2k -i R 2g+1 (Z) 

The differentials fulfill 

(3-7) / Cj=Sk,j, / Cj=Tk,j, Tk,j=Tj,k, j,k = l,...,g. 
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Let us now pick g numbers (the Dirichlet eigenvalues) 
(3-8) (£i)?=i = G^°-i)?=i 

whose projections lie in the spectral gaps, that is, /ij £ [E^-i, E 2 j], j — 1, ■ ■ • , g- 
Associated with these numbers is the divisor 

and we can define g numbers {pt>j{x,t)) a - 1 — (fij(x,t),<Tj(x ) t))j_ 1 via Jacobi's 
inversion theorem by setting 

(3.10) « £o (%,,t)) = fi^A) + J^o + 12^ 2 

such that /ij(Q,0) = /X/. Here U and C/ 2 denote the 6-periods of the Abelian 
differentials w Poc ,o an d u Poo ,2, respectively, defined below, and A Eq (a Eg ) is Abel's 
map (for divisors). The hat indicates that we regard it as a (single- valued) map 
from Kg (the fundamental polygon associated with K g by cutting along the a and 
b cycles) to C 9 . 

Next we introduce 

(3.11) z(p,x,t) = E Eo -A Eo (p) + a Eo {V k(Xtt) ) e C 9 , 
where S Bo is the vector of Riemann constants 

/oio\ ~ _ i + z)fc=i T j,fc 

(3-12) "E ,j= 2 ' J — 

By we denote the Riemann theta function associated with /C g defined by 

(3.13) 6(z) = gx p 2 ^ + ^-'^-M , zeC 9 . 

Note that the function 0(z_(p, x,t)) has precisely g zeros £ij(x,t). This follows from 
Riemann's vanishing theorem (cf. [33J Theorem A. 13]). 

Introduce the time-dependent Baker- Akhiezer function 
(3.14) 

x ' * } = e(z( P ^x,t)) e(z(p,o,o)) exp v " ffl -L ^° - 12lt L " P -V ■ 



En J Ei 



Here w Pooi o & n d w Pooi 2 are normalized Abelian differentials of the second kind with 
a single pole at p^, and principal part w~ 2 dw and w~ i dw in the chart w(p) = 
iiz^ 1 / 2 for p = (z,±), respectively. The Abelian differentials are normalized to 
have vanishing aj periods and have the following expressions 



(3-15) u Poot0 = ^ -j/2 dir, 



l n^i^- A,) 
2i i?y 2 . 



with Xj € (E 2 j-i,E 2j ), j = l,...,g, and 

(3-16) w Poo , 2 = -jya dir, 

H 2g+1 
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where Xj,j = 0, . . . , g, have to be chosen such that they fulfill J2j=o A? = 1 j=o ■ 
We also remark 

(3.17) iP q (p,x,t)i> q (p*,x,t) = T[ — ^M , p=(z,±). 

i : ''' 

Then our background KdV solution is given by 

a 

(3.18) V q {x,t) = E + J2(E2j-i+E 2j - 2A,-) - 29^ ln0(z( Poo , x, i)) . 

3=1 

The Abelian differentials of the third kind u qi q2 , with simple poles at q± and 
(j2, corresponding residues +1 and —1, vanishing a-periods, and holomorphic on 
K g \ {qi, q 2 }, are explicitly given by (;16, Appendix B]) 

V 2(7r-7r(pi)) 2(7r-7r(p 2 )j ^+1 

(3 - 20) 2(,-n( Pl )) + PlP ~ {Z) JR%r' 

where pi,p 2 € ^3 \ {Poo} and -P PlP2 (z) and P PlPoo (z) are polynomials of degree 
g — 1 which have to be determined from the normalization J a w Pl P2 = and 
J a uj Pi Poo = 0, respectively. In particular, 

(3.21) u, pp , = (_^ + P w » (7r )j 
We will also need the Blaschke factor 

(3.22) B(p,p) = exp f / w p p « ) = exp ( f u pp *\ tt(p) € E, 

where -E(p) is _E if p < -Bo, either i?2j-i or i? 2 j if p G (E 2 j-i,E 2 j), 1 < j < g. It 
is a multivalued function with a simple zero at p and simple pole at p* satisfying 
\B(p,p)\ = 1, p € 9II + . It is real- valued for 7r(p) € (— 00, Eo) and satisfies 

(3.23) B(E ,p) = 1 and B(p*, p) = B(p, p*) = B(p, p)- 1 
(see e.g., [35]). 

The Baker- Akhiezer function is a meromorphic function on JC g \ {poo} with an 
essential singularity at p^. The two branches are denoted by 

(3.24) tp q>± (z,x,t) = i/j q (p,x,t), p=(z,±), 
and it satisfies 



H q (t)ip q (p,x,t) = Tr(p)ij> q (p,x,t), 
d_ 
dt' 



(3.25) -i> q (p,x,t) = P q , 2 (t)^ q (p,x,t). 
Here 

H q (t) = d 2 x + V q (.,t), 

(3.26) P q , 2 (t) = -id x + 6V q (.,t)d x + 3V q , x (.,t), 
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are the operators from the Lax pair for the KdV equation, that is, 

(3-27) jH q {t) = H q {t)P q , 2 {t) - P q , 2 {t)H q {t). 

It is well known that the spectrum of H q (t) is time independent and consists of 
g + 1 bands 

g-l 

(3.28) a(H q {t)) = |J [E 2j ,E 2j+1 ] U [E 2g , oo). 

3=0 

For further information and proofs we refer to |16j . 

4. The Inverse scattering transform and the Riemann-Hilbert 

PROBLEM 

In this section we recall some basic facts from the inverse scattering transform 
for our setting. For further background and proofs we refer to [4], [14], and [12] 
(see also [25]). 

Let tp qt ±(z, x, t) be the branches of the Baker- Akhiezer function defined in the 
previous section. Let ^±(z,x,t) be the Jost functions for the perturbed problem 

(4.1) (-d£ + V(x,t))i/>±(z,x,t)=zil) ± (z,x,t), 
defined by the asymptotic normalization 

(4.2) lim e^ fc W ty ± (z, x, t) - W± {z, x, t)) = 0, 
where k(z) denotes the quasimomentum map 

(4.3) k(z) = - [ w Po0! o, p=(z,+). 

J Eq 

The asymptotics of the two projections of the Jost function are (cf. 2U Theo- 
rem 2.3]) 

/ f ±oa 1 \ 

(4.4) iP±(z,x,t) = rp qt± (z,x,t)(lT J (V -V q )(y,t)dy^= + o{l/^)j, 

as z — > oo. We will assume that the poles of the Baker- Akhiezer function ^ are 
all different from the eigenvalues pj without loss of generality (otherwise just shift 
the base point (x ,t Q ) — (0,0)). 
One has the scattering relations 



(4.5) T(z)tp T (z,x,t) = ip±(z,x,t) + R±(z)ip±(z,x,t), z £ a(H q ), 

where T(z), R±(z) are the transmission respectively reflection coefficients. Here 
ip±(z,x,t) is defined such that ip±(z,x,t) — lim £ iQip±(z + ie,x,t), z e a(H q ). If 
we take the limit from the other side we have ip± (z, x, t) = lim e ^o "0± ( z — i £ j x i t)- 

The transmission and reflection coefficients have the following well-known prop- 
erties: 

Lemma 4.1. The transmission coefficient T(z) has a meromorphic extension to 
C\a(H q ) with simple poles at the eigenvalues pj. The residues ofT(z) are given by 



(4.6) Res Pj T(z) 



2-^23+1 (Pj) l±.j 

ftlw "IEU(pj-w)^ 



3 
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wh 



ere 



(4-7) 7ij= / \1>±(j>s,V,t)fdy 

are referred to as norming constants and ip-(pj,x,t) = Cjip+(pj,x,t). 
Moreover, 



(4.8) T(z)R+(z)+T(z)R-{z)=0, \T(z)\ 2 + \R ± (z)\ 2 = 1. 

In particular, one reflection coefficient, say R{z) — R + {z), and one set of norming 
constants, say jj = 7+j , will be sufficient for us. 

We will define a sectionally meromorphic vector on the Riemann surface 1C g as 
follows: 

!(T(~\ TP-( Z ' X ' 1 ) ji+(z,x,t) \ ( J 

) I I *\ If U \ 

KjZ^T) 1 y z h g ,-(z,x,t)) > p -(*>-; 

We are interested in the jump condition of m(p,x,t) on E, the boundary of fl± 
(oriented counterclockwise when viewed from top sheet 11+ ). It consists of two 
copies E± of cr(H q ) which correspond to non-tangential limits from p — (z, +) with 
±Im(z) > 0, respectively to non-tangential limits fromp = (z, — ) with =plm(z) > 0. 

To formulate our jump condition we use the following convention: When repre- 
senting functions on E, the lower subscript denotes the non-tangential limit from 
fl + or H_, respectively, 

(4.10) rn±(po) = lim m(p), p € E. 

n±9p->j>o 

Using the notation above implicitly assumes that these limits exist in the sense that 
m(p) extends to a continuous function on the boundary away from the band edges. 
Moreover, we will also use symmetries with respect to the sheet exchange map 



(4.11) p* 
and complex conjugation 

(4.12) p = 



I(z,=f) fcrp = (z,±), 
[poo for p = Poo, 

(z,±) forp=(z,±)£E, 
(z,T) forp= (z,±) G E, 
Poo for p = Poo- 



In particular, we have p = p* for p G E. 

Note that we have m±(p) — m^(p*) for rh(p) = m(p*) (since * reverses the 
orientation of E) and m±(p) = m±(p*) for rh(p) = m(p). 

Note that we have the following asymptotic behavior for m(p, x, t) near poo- 

1 f°° 1 
(4.13) m(p) = (1 l)-x^J x (V-V q )(y,t)dy(-1 l)+o(-^), P=(z,±) 

for p near poo ■ Here we made use of (|4.4[) and 



If 00 1 
(4-14) T{z) = l + — T {V - Vq){ y 7 t)dy + o(^) 

(cf. [2H Corollary 3.7]). 

We are now ready to derive the main vector Riemann-Hilbert problem: 
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Theorem 4.2 (Vector Riemann-Hilbert problem). Let S+{H(0)) = {R(X), A £ 
o-(H q ); (pj, 7j), 1 < j < JV} the right scattering data of the operator H(0). Then 
m(p) = m(p, x, t) defined in (|4.9[) is meromorphic away from S and satisfies: 

(i) The jump condition 

(4-15) ' 

for 

(ii) the divisor 

(4.16) (mi) > -!>/}( x , t ). - P £ , (m 2 ) > -£>/i(x,t) ~ ^p* 

and poZe conditions 

(4.17) 

,1/2 



/ , x , ~ 2 - R 2g+l(ft) 7j Tpq(P,X,t) \ 

v mib) + nLi(^ - ^ - * ^w) m2(p) J - "^ ( -* } " nearpj ' 

1 /2 

Irfff — ? r? IxT m iW + m 2W J >-^Kx,t), near p^ 

(iii) i/ie symmetry condition 

(4.18) m(p*)=m(p)(j J) 

(iv) and the normalization 

(4.19) m( Poo ) = (1 1) . 
Here (/) denotes the divisor of f and 

(4.20) 2? £ = 5>«' %=E^" 

denotes the divisor corresponding to the points pj = [pj,+) £ IC g . The phase <fi is 
given by 

(4.21) 0(p, -) = -24i / cj Poo , 2 -2i- w Poc , e« /orpeE. 

1 •'Po 1 •'Po 

Moreover, we have set 



(4 22) e(pj;f) _ gU(P,»,*))gU(P*,0,0)) 
such that 

Here we have extended our definition of R to E such that it is equal to R(z) 
on E_|_ and equal to R(z) on E_. In particular, the condition on E+ is just the 
complex conjugate of the one on E_ since we have R(p*) — R{p) and m± (p*,x,t) = 
m±(p, x, t) for p £ E. 
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Proof. The jump condition follows by using (|4.5|) and ()4.8[) . By Riemann's vanish- 
ing theorem (cf. [HH Theorem A. 13]) the Baker- Akhiczcr function ip q has simple 
zeros at jlj(x,t) and simple poles at p,j, j = 1, ... ,g. Moreover, the transmission 
coefficient T(z) has simple poles at the eigenvalues pj, j = 1,...,N. Thus the 
divisor conditions (|4. 16[) are indeed fulfilled. The pole conditions follow from the 
fact that the transmission coefficient T[z) is meromorphic in C\ o~(H q ) with simple 
poles at pj and its residues are given by (|4.6j) . The symmetry condition (|4.18p 
obviously holds by the definition of the function m(p). The normalization (|4.f 9[) is 
immediately clear from (|4.13p . □ 



We note that the symmetry condition is in fact crucial to guarantee that the 
solution of this vector Riemann-Hilbert problem is unique. 

Theorem 4.3. The vector m(p) defined in (14. 9|) is the only solution of the vector 
Riemann-Hilbert problem (|4.15|) - (l4.19p . 

Proof. The argument is similar to [22j Thm. B.l]. It suffices to show that the corre- 
sponding vanishing Riemann-Hilbert problem, where the normalization condition 
(|4.19[) is replaced by m(poo) = (0 0), has only the trivial solution. 

Let rh be some solution of the vanishing Riemann-Hilbert problem. We want 
to apply Cauchy's integral theorem to ih(p)fw (jf). To handle the poles of m 
we will multiply it by a meromorphic differential c?f2 which has zeros at p and 
p* and a simple pole at p^ such that finally the differential rh{p)fh) (p*)dQ(p) is 
holomorphic away from the contour. Here rh^ denotes the adjoint (transpose and 
complex conjugate) vector of rh. 

More precisely, let 

(4.23) <m = nf=ifr-Mi) d7r 

pV^ 

and note that — (Ylj(z— Pj)) R2g+i( z ) is a Herglotz function. That is, it has positive 
imaginary part in the upper half-plane (and it is purely imaginary on o~(H q )). Hence 
rh{p)rh T (p)dQ(p) will be positive on S. 
Next, consider the integral 



(4.24) 0= / m(p)m<(F)dn(p), 



D 

where D is a ^-invariant contour consisting of two loops on the upper and on the 
lower sheet encircling none of the poles pj , p* . We first deform D to a *- invariant 
contour consisting of several parts: Two pieces D± wrapping around the ± side of 
S plus a number of small circles -D+j, D-.j around the poles pj, p*, respectively. 
Then the contribution from E is given by 



m(p)ml (p*)d£l(p) = / (m + (p)m ! _(p*) + ■m^(p)m\_(p*))dfl(p) 
(4.25) = / m-(p){J(p) + J f (p*))mL(p*)dfl(p) > 
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and the contribution from the poles is given by 

m (p) rfv (p* ) dfl (p) 



uf =1 (^+, 3 uD_, 3 ) 

N 

= ^ ^Res Pj . rh(p)m\p*)dS}(p) + Res p * rh(p)m\p*)drt(p)j 

N 

(4.26) = 2^ Resp. m{p)m\p*)dtt(p). 

3=1 

To compute the residues we use the pole conditions (|4.17p which imply (using 

©33)) 



Res Pj m{p)m\p*)dn{p) = ^ - a ^ q { Pj ) 2 m 2 { Pj f > 0. 



In particular, both contributions to the integral (I4.24p are non-negative and thus 
both must vanish. It follows from the that rh = vanishes along £ and consequently 
rh(p) — as desired. □ 

We will also need another asymptotic relation 
(4.27) mi • m 2 = 1 + (V - V q )(x, t)^ + o^" 1 ). 

which is immediate from the following well-known result. 
Lemma 4.4. We have 

(428) T(z) ^-(,,x,0^ + (^x,t) =1+ 2^-^^^z + ° (z } 

Proof. We will use the following representation of the Jost solutions 

(4.29) ip±(z,x,t) = il) q> ±(z,x,t)exp ( =p / (m±(z,y,t) - m q ,±(z,y,t))dy), 
where 

( ±\ / j \ ,^,± 

m±(z,a;,<) = ±— — , m qt± {z,x,t) = ±-^— 

ip±(z,x,t) ip q ,±(z,x,t) 

are the Weyl-Titchmarsh functions. Here the prime denotes differentiation with 

respect to x. Using the expansion of the Weyl m-functions (cf. [29l Lemma 6.1]) 

and the one for logT(z) (cf. [29l Theorem 6.2]) for z — > oo proves the claim. □ 

For our further analysis it will be convenient to rewrite the pole conditions as 
jump conditions following the idea of Deift, Kamvissis, Kriecherbauer, and Zhou 
|11) . For that purpose we choose e so small that the discs \ir(p) — pj\ < e are inside 
the upper sheet II + and do not intersect with the spectral bands. Then redefine 
m(p) in a neighborhood of pj respectively p* in the following way: 

™>{P) I l3 ( P ,x,t) l I , pen+ > 

\ ir(p)-pj / 
/ 1 7j(p ,x,t) \ 

m(p){ ' '• " . 



(4.30) m(p) = < 



pen 

^m(p), else, 
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where Jj(p, x, t) is a function which is analytic in < \tt(p) — pj\ < £, p G 11+ and 
satisfies 

lim T (p xt) ^^± = ggfc^ T 
For example, we can choose 



Lemma 4.5. Suppose m{p) is redefined as in (I4.30p . TTien m(p) is meromor- 
phic away from E anc? satisfies (|4.15l) . (|4.18l) . (|4.19[) . £/ie divisor conditions change 
according to 

(4-31) (mi) > -Vypjy, (m 2 ) > -T>fi{x,t) 

and the pole conditions are replaced by the jump conditions 



(4.32) 



/ 1 

m+(p) = m_(p) 7 j(P,^,«) x 

\T(P)-Pj 

m+(p) = m-(p) I t(p)-pj 




(4.33) E s (p) ={?en±: Ivr(g) - z| = e}, p = (z, ±), 

is a small circle around p on the same sheet as p. It is oriented counterclockwise 
on the upper sheet and clockwise on the lower sheet. 

Proof. Everything except for the pole conditions follows as in the proof of Theo- 
rem l4.2l That the pole conditions (|4.17j) are indeed replaced by the jump conditions 
(|4.32[) as m (p) is redefined as in (|4.30j) can be shown by a straightforward calcula- 
tion. □ 

The next thing we will do will be to deduce the one-soliton solution of our 
Riemann-Hilbert problem, i.e., the solution in the case where only one eigenvalue 
p corresponding to one bound state is present and the reflection coefficient R(p) 
vanishes identically on K g . 

Lemma 4.6 (Onc-soliton solution). Suppose there is only one eigenvalue and a 
vanishing reflection coefficient, that is, S + {H{t)) = {R(p) = 0, p G E; (p, 7)}. Let 



(4.34) c q ,y(p,x,t) = l + "/W( x> t)(if> q {p,x,t),il; q (p,x,t)) = I + 7 / ip q (p,y,t) dy 

J X 

and 

(4.35) g7feg ,t) = ^(p >8>t ) + ^^^ 

z-p C q .j{p,X,t) 

Here the dot denotes a derivate with respect to p and Wi x ^(f, g) = (f(x)g'(x) — 
f'(x)g(x)) is the usual Wronski determinant, where the prime denotes the derivative 
with respect to x. 
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(4.37) (V - V q )(x, t) = log (c qn (p, x, t)) . 



Then the unique solution of the Riemann-Hilbert problem (|4.15j) — (|4. 19[) is given 

by 

mo(p) = (/(p*,x,t) f(p,x,t)), f{ P , x ,t) = ^^ X f 

lj} q (P,X,t) 

In particular 

(4.36) J (V -V q )(y,t)dy = -2— log (c qn (p,x,t)), 

or 

d* 
dx 2 

Proof. Since we assume the reflection coefficient vanishes, the jump along E disap- 
pears. Moreover, since the symmetry condition (|4.18p has to be satisfied it follows 
that the solution of the Riemann-Hilbert problem (|4.15[) - (14.19l) has to be of the 
form ma(p) = (f(p* , x, t) f(p, x, £)) . The divisor conditions (|4.16j) follow from the 
fact that the Baker- Akhiezer function ip q has simple zeros at fij(x,t) and simple 
poles at fij, j = l,...,g and by construction of ip q .-y . It is obvious that the nor- 
malization condition (|4.19l) holds. Thus it is only left to check the pole conditions 
(|4.17p . For this purpose we compute 

Hm(z- p)f{p*) = ' ■ W( x ^(ijj q (p,x,t),'ip q (p* ,x,t)) 

p—>-p c q ^(p, , X, t) 

j(p,x,t) 2R^ 2 +1 (p) 

c qn (p,,x,t) n*=i(p-Mfc)' 

where we defined 

7 (p, x, t) = 7 y^4 = x > f ) e * 0(P) > 

IpqiP ,X,t) 

and we used (cf. [H Equ. (1.87)]) 

2i? 1/2 (z) 
WW q , T (zU q , ± (z))=± Tig 29+11 



T[fc=i( z -Mfc) 



Moreover, 



,. ,/ v , , 7 ,. W( x>t )(ip q (p,x,t),il> q (p,x,t)) 
hm f{p) = 1H hm 

P^P C q ^{p,X,t) Z ~P 

1 7 V , , ,\ y ^ q (P,X,t)-ll}'(p,X,t) 

= 1 H 7 77 Vw, a;, *) lim — 

c qn {p,x,t) L p^p z — p 

- ip Jp, x, t) hm -2 2 

* p^p z — p J 

= 1 + 7 7^W( x>t )(lp q (p,X,t), tp q (p,X,t)) 



Hence we see that the pole conditions (|4.17j) are satisfied. 

The formula (|4.36j) follows after expanding around p — p^, that is, 

f(p,x,t) = 1 + , - , 7 , —ip q (p,x,t)(i(j q (p 7 x 7 t)m q (p 7 x 7 t) - ip' q (p,x,t)j 

(z p)c q ^{p,x,t) \ / 
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where we have used that the Weyl-Titchmarsh m-function has the following as- 
ymptotic expansion for p near poc (cf. [29[ Lemma 6.1]) 

(4.38) m q , ± (z,x,t) = = ±iy r z+ Y^L + (z~% P =(z,±). 

Thus comparing with (|4.13p proves the equation (|4.36[) . 

To see uniqueness, let rho(p) be a second solution which must be of the form 
fho(p) — {f{p*) f(p)) by the symmetry condition. Since the divisor T)p,( x ,t) is 
nonspecial, the Riemann-Roch theorem implies f(p) = af(p)+/3 for some a, ft € C. 
But the pole condition implies /J = and the normalization condition implies 
a = 1. □ 

Since up to quasi-periodic factors ip q (p,x,t) is a function of x — v{p)t, where 

-12Ref^ + V, 2 



(4.39) v(p) 



Re /En' +)a; Poo,0 



we will call v(p) the velocity of the corresponding soliton. 

5. The stationary phase points and the nonlinear dispersion relation 

In this section we want to look at the relation between the energy A of the 
underlying Lax operator H q and the propagation speed at which the corresponding 
parts of KdV solutions travel, that is, the analog of the classical dispersion relation. 
If we set 

-12Re(i J^ +i£ ' +) ^ a ) -12Im/j, A+iEl+) a, Poo , 2 

(5.1) w(A) = lim \H° ^llZ = lim i^o 

e->0 D./; r(A+ie,+) \ e ^ t f (A+is,+) 

our first aim is to show that the nonlinear dispersion relation is given by 

(5.2) V ^ = T 

Recall that the Abelian differentials are given by ()3.15[) and (|3.16p . 

For p 6 M\a(H q ) the denominator is nonzero and the formula agrees with the 
soliton velocity defined in (|4.39|) . In particular, recalling the definition of our phase 
<j> from (|4.2ip . this implies 

(5.3) = 7 IW(p,~) = 

in this case. In particular, this definition reduces precisely to the definition of the 
velocity of a soliton corresponding to the eigenvalue p (cf. the discussion after 
Lemma |4.6[) . 

For A € cr(H q ) both numerator and denominator vanish on o{H q ) by (|3.2p . 
Hence by virtue of the de l'Hospital rule we get 

(5 ' 4) VW - - II^A-A,) ■ 

that is, 

(5.5) u ( A ) = f ^ 0'(A,f) = O. 

In other words, v(X) coincides with a stationary phase point in this case. 
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So let us discuss the stationary phase points, that is the solutions of </>'(A, § ) = 0, 
next. The solutions are given by the zeros of the polynomial 

(5.6) nfliz-x^ + ^fliz-x,). 

3=0 j=l 

Since our Abelian differentials are all normalized to have vanishing a., -periods, the 
numbers Xj, < j < <?, are real and different with precisely one lying in each 
spectral gap, say A.,- in the j'th gap. Similarly, Xj, < j < g, are real and different 
and Xj, 1 < j < <?, sits in the j'th gap. However Ao can be anywhere (see [34j 
Sect. 13.5]). 

The following lemma clarifies the dependence of the stationary phase points on 
x/t. 

Lemma 5.1. Denote by Zj(v), < j < g, the stationary phase points, where 
v = x/t. Set Xq = — oo and A g +i = oo, then 

(5.7) Xj < Zj[v) < Xj+i 

and there is always at least one stationary phase point in each spectral gap. More- 
over, Zj(v) is strictly decreasing with 

(5.8) lim Zj(v) = Aj+i and lim Zj(v) = Xj. 



-oo 



Proof. Since the Abelian differential ^> Poo .2 + vuj Pix .q has vanishing a periods, the 
polynomial (|5.6[) must change sign in each gap except the lowest. Consequently 
there is at least one stationary phase point in each gap except the lowest, and they 
are all different. Furthermore, by the implicit function theorem, 

/ = g(*j) = IlLi( z j ~ A fc) 

q'(zj) + vq'{zj) 12 flLo,^ ( z i ~ **) ' 

where 

g g 

q{z) = 12 J[{z-X k ), q(z) = l[(z-X k ). 
k=a fc=i 
Since the points Afe are fixed points of this ordinary first order differential equation 
(note that the denominator cannot vanish since the Zj's are always different), the 
numbers Zj cannot cross these points. Combining the behavior as v — > ±oo with 
the fact that there must always be at least one of them in each gap, we conclude 
that Zj must stay between Xj and Aj+i. This also shows z'j < and thus Zj(v) is 
strictly decreasing. □ 



In other words Lemma 15.11 tells us the following: As v = x/t runs from — oo 
to oo we start with z g (v) coming from oo towards Ei g , while the other stationary 
phase points Zj, j = 0, . . . ,g — 1, stay in their spectral gaps until z g (v) has passed 
Ei g and therefore left the first spectral band \E<z g ,oo). After this has happened, 
the next stationary phase point z g -\(v) can leave its gap [Ei g —\,Ei^) while z g (v) 
remains there, traverses the next spectral band [I?2 fl -2> ^g-i] and so on. Finally 
z (v) traverses the last spectral band [E ,Ei] and moves to — oo. So, depending 
on x/t there is at most one single stationary phase point belonging to the union of 
the bands cr(H q ) which is then the one solving (|5.5[) . 

Lemma 5.2. The junction v(X) defined in (|5.ip is continuous and strictly decreas- 
ing. Moreover, it is a bijection from R to R. 
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Proof. That v(X) defined in (|5.1|) is continuous is obvious except at the band edges 
A = Ej. However, in this case (|5.1|) becomes (|5.4|) by using the de l'Hospital rule. 
The function v(X) defined in (|5.4[) is obviously continuous at the band edges Ej 
since Xj lies in the j'th gap and thus does not hit the band edges. 
Furthermore, for large A we have 

(5.9) hm ^ = 1, hm = 1, 

V ' A^-oo -4A A^+oo -12A 

which shows liniA->±oo v(X) = =pco. 

In the regions where there is one stationary phase point Zj(v) G a(H q ) we know 
that Zj(v) is the inverse of v(X) and monotonicity follows from Lemma 15. II In the 
other regions we compute 



(5.10) v'(z) 



n- D 1 / 2 / \ f ( z >+) 



Let the stationary phase points be ordered such that we have Zj(v(z)) = z for 
z G {Ezj—i — e,E2j—i] and Zj—i(v(z)) = z for z G [E2j,E2j + e). Then we claim 
Zj(v(z)) < z < Zj-x(v(z)) for z G (E2j-i, E%j) (set Z-\ — oo). In fact, the above 
differential equation implies that v(z) can cross the curve Zj{v{zj) only from below 
and hence must stay above this curve since it starts on this curve at z — £<y_i. 
Similarly, it can cross the curve Zj-i(v(z)) only from below and would hence remain 
above this curve afterwards. Thus this can only happen at z = E2j. □ 

In summary, we can define a function £(x/t) via 

(5.11) v(C) = j. 

In particular, different solitons travel at different speeds and don't collide with each 
other or the parts corresponding to the continuous spectrum. 

Moreover, there is some Co f° r which w(Co) — and hence there can be stationary 
solitons provided Co ^ °~{H q ). 

Lemma 5.3 (Stationary solitons). There exists a unique Co such that u(Co) = 0. 
Moreover, if Co € cr(H q ) or Xq 6 cr(H q ), then Co = Ao- In particular, Co S °~(Hq) if 
and only if Xq G a(H q ). 

Proof. Existence and uniqueness of Co follows since v is a bijection. It is left to 
show that Co = Ao if C G a {H q ) or Ao € a{H q ). Assume Co G ^{Hq)- Then using 
u(Co) = and pT3]l we get 

9 9 

IICCo - Aj) = (Co - A ) n(Co - Xj) = 0. 

J=0 i=l 

Since Aj G (E2j-i,E2j), j = l,...,g it follows Co = Ao- Now suppose A G o-(H q ) 
and again use (I5.4[) to get 

KAo)I[(Ao-A i ) = 0. 

i=i 

Since Aj G (E^j-i, ^2j)j J ' = 1' • ■ • ; 9 we obtain w(Ao) = and thus Co = Ao- □ 
In summary we conclude that depending on v = | there can occur three cases: 
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Figure 1. The lens contour near a band containing a stationary 
phase point Zj and its flipping image containing z*. Views from 
the top and bottom sheet. Dotted curves lie in the bottom sheet. 



(i) ((v) £ (E2j,E 2 j+i) for some j = 0,...,g (setting £23+1 = oo). In this 
case £(v) is a stationary phase point and all other stationary phase points 
lie in open gaps. 

(ii) £(u) € R \ cr(H q ) and all stationary phase points lie in open gaps. 

(iii) £(u) = Ej f° r some j and all other stationary phase points lie in open 
gaps. 

These three cases define corresponding regions in the (a;, t)-plane: the oscillatory 
region (case (i)), the soliton region (case (ii)), and the transitional region (case (iii)). 

Case (i): The oscillatory region. Note that in this case we have 
(5.12) cj) { Zj )/i = 1/2 ' ; > 0. 

lR 2g+l( Z j) 

Suppose C(v) — Zj(v), belongs to the interior of the band [E^j, E%j+i] (with E 2g+ i = 
oo). Then we introduce the "lens" contour near that band as shown in Figure [TJ 

The oriented paths Cj = Cji U Cj2, C* — C* x U C| 2 are meant to be close to the 
band [E 2j ,E2 J+1 }. 

Concerning the other bands [-E^fej -E/2fe+i]; k ^ j, k = 0, . . . , g (setting E 2g+ i = 
oo), one simply constructs "lens" contours near each of the other bands [^fe, -E^fc+i] 
and [E% k , E^ k+1 ] as shown in Figure [2J 

The oriented paths Ck,C k are meant to be close to the band [E 2 k, E 2 k+i]- In 
particular, these loops must not contain any of the eigenvalues pj. 

Then an investigation of the sign of TLe(4>) shows the following: 

|Re(0(p)) > 0, pe Djt U D k , ir(p) < £(x/t), 
\Re(<f>(p)) < 0, pe D j2 U D k , n(p) > Q(x/t), 

with k = l,...,g, k^j. 
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E2k 









Figure 2. The lens contour near a band not including any sta- 
tionary phase point. Views from the top and bottom sheet. 

Observe that our original jump matrix (|4.15|) has the following important fac- 
torization 

(5.13) J(p)^b^(p)- 1 b + (p), 
where 

(5.14) b-(p)={ Q > \ > j, b + (p) = [ R(j))Q{p)etm 

which is the right factorization for p g T,\C(x/t) — £ n 7r _1 ((C(x/t), oo)), i.e., 
ir(p) > Q{x/t). Similarly, we have 

(5.15) J {p )=B_ {p )^( 1 -\f^ _±_) B+{P ), 

V U i-\R(p)\ 2 J 



where 



1 0\ /i fl(p*)e(p-)c- t ^'' 



(5.16) B-(p) = R{p) e( P )c^ , , B+(p) = x ^™ 

V i-ii!wi» V V° 1 / 

This constitutes the right factorization for p £ C(x/t) = Sn7r _1 ((— oo, ((x/t j), i.e., 
w(p) < ((x/t). Here we have used R(p) = R(p*), for peS. To get rid of the diag- 
onal part in the factorization corresponding to ir(p) < ((x/t) and to conjugate the 
jumps near the eigenvalues we need to find the solution of the corresponding scalar 
Riemann-Hilbert problem, the so-called partial transmission coefficient. Again we 
seek a meromorphic solution. This means that the poles of the scalar Riemann- 
Hilbert problem will be added to the resulting Riemann-Hilbert problem. On the 
other hand, a pole structure similar to the one of m is crucial for uniqueness. We 
will address this problem by choosing the poles of the scalar problem in such a way 
that its zeros cancel the poles of m. The right choice will turn out to be T>o_ (that 
is, the Dirichlet divisor corresponding to the limiting lattice defined in (|2.8|l ). 
Define a divisor T>p/ X ^ of degree g via 

(5-17) Oso =°Le ( t> AOM)) +£0/*)) 
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where 

(5.18) S e (x/t) = -2 A E( Pk )M) + i^l log(l - |i?| 2 )0, 

with C{x/t) = Sfl 7r _1 ((-oo,C(x/*)) and ((x/t) as defined in ([5~TT]) . 

Then "D^u.^ is nonspecial and Tr(£ , j(x,t)) = Vj(x,t) G K with precisely one in 
each spectral gap (see [22]). 

We define the partial transmission coefficient as 
(5.19) 

_ g(z(p 00 ,x,t) + £(x/t)) 6>(z(p,x,t)) 

0U(Poo,x,t)) 0(«(p, a, *)+£(*/*))' 

where <5(x, t) is defined in (j5 . 18[) and Cc^ P2 is the Abelian differential of the third 
kind with poles at p\ and p2- 

The function T(p,x,t) is meromorphic in K, g \ C(x/t) with first order poles at 
Pk < C(x/t), i>j(x,t) and first order zeros at (Mj(x, t). 

Lemma 5.4. T(p,x,t) satisfies the following scalar meromorphic Riemann-Hilbert 
problem: 

T + (p,x,t) = T_(p,x,i)(l - |i?(p)| 2 ), p G C(x/t), 



(5.20) {T(p,x,t))= 2_j V pI~ P ft+%M)-%M> 

p h <C(x/t) Pk <c(x/t) 

T(poa,X,t) = 1. 



Moreover, 
(i) 

T(p*, x, tJTfo x, i) = f[ Z -^ X '^ , z = n(p). 

7=1 Z ^A 25 ' f J 



(ii) T(p, x, t) = T(j>, x, i) and in particular T{p, x, t) is real-valued for tt(j>) G 
R\a(H q ). 

Proof. The argument is similar to [HI Thm. 4.3]. The solution of a Riemann- 
Hilbert problem on the Riemann sphere is given by the Plemelj-Sokhotsky formula. 
Since our problem is now set on the Riemann surface JC g the Cauchy kernel is given 
by the Abelian differential of the third kind w pPoc (cf. [35]). In particular, T(p, x, t) 
satisfies the jump condition from (|5.20p along C(x/t). Next we have to check that 
the function T(p, x, t) extends to a single- valued function on JC g . For that purpose 
note that the only possible contribution which causes multi-valuedness may come 
from the 6-cycles since all Abelian differentials are normalized to have vanishing 
a-periods. So for the 6^-periods £= 1, ... ,g we compute for p G C(x/t) 



e±oT(p-ie,x,t) V J c(x/t) pk ^ x/t) I 

which is indeed equals 1 by the choice of Se in (|5.18[) . 
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Concerning the poles and zeros of the function T(j>, x, t) we see that by Riemann's 
vanishing theorem (cf. 34, Theorem A. 13]) and the choice of the divisor T>pt Xtt ) 
defined by (|5.17[) that the ratio of theta functions is meromorphic with simple zeros 
at pLj and simple poles at Moreover, from the product of the Blaschke factors 
we get that T has simple poles at pk and simple zeros at p* k for which p^ < £(x/i) 
is valid. 

To prove uniqueness let T be a second solution and consider T/T. Then T/T 
has no jump and the Schwarz reflection principle implies that it extends to a mero- 
morphic function on K, g . Since the poles of T cancel the poles of T, its divisor 
satisfies (T/T) > —T>^ x t y Since "Dp,f x ,t) i s nonspecial, T/T has to be a constant 
by the Riemann-Roch theorem (cf. [3~U Theorem A.2]). Setting p = p^, we see 
that this constant is one, that is, T = T as claimed. 

Finally, T(p,x,t) = T(j5,x,t) follows from uniqueness since both functions solve 

We will also need the expansion around p^ given by 

Lemma 5.5. The asymptotic expansion of the partial transmission coefficient for 
p near poo is given by 

(5.21) T(p,x,t) = l±^^-+0(h p=(z,±), 



z 



where 



(5.22) 



T 1 (x,t)=- V 2/ u; Pao . + — log(l - |iJ| a )«p„ 



id y In 




Q{z(poo,y,t)) 

and w Poo ,o * s ^ e Abelian differential of the second kind defined in (|3.15|) . 

Proof. This can be verified similarly as in the case of the full transmission coefficient 
(cf. [221 Theorems 6.2 and 6.3] and by expanding the ratio of theta functions near 

Poo- □ 

Now that we have solved the scalar Riemann-Hilbert problem for T(p,x,t) we 
can conjugate our original Riemann-Hilbert problem. 

Since to each discrete eigenvalue there corresponds a soliton, it follows that 
solitons are represented in our Riemann-Hilbert problem by the pole conditions 
(|4.32jl . For this reason we will study how poles can be dealt with in this section. 
We will follow closely the presentation of [24] Section 4] . 

In order to remove the poles there are two cases to distinguish. If pj > £(x/t), 
the jump at pj is exponentially close to the identity and there is nothing to do. 

Otherwise, if pj < Q(x/t), we need to use conjugation to turn the jumps at these 
poles into exponentially decaying ones, following [TT]. It turns out that we will have 
to handle the poles at pj and p* in one step in order to preserve symmetry and in 
order to not add additional poles elsewhere. 

Moreover, the conjugation of the Riemann-Hilbert problem also serves another 
purpose, namely that the jump matrix can be separated into two matrices, one 
containing an off-diagonal term with exp(— tip) and the other with exp(f0). Without 
conjugation this is not possible for the jump on C(x/f) — E n 7r _1 ((— oo, ((x/t))), 
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since in this case there also appears a diagonal matrix if one wants to separate the 
jump matrix. 

For easy reference we note the following result, which can be verified by a 
straightforward calculation. 



Lemma 5.6 (Conjugation). Assume that ECS. Let D be a matrix of the form 
(5.23) D(p) 



d(p*) 
dip)) ' 

where d : K, g \S — > C is a sectionally analytic function. Set 

(5.24) fhip) = m(p)Dip), 
then the jump matrix transforms according to 

(5.25) J{p) = D-{p)- 1 Jip)D+ip). 

rh(p) satisfies the symmetry condition (|4.18p if and only if m(p) does. Further- 
more, rhip) satisfies the normalization condition (|4.19l) if m(p) satisfies (|4.19l) and 
d(Poo) = 1- 

Lemma 5.7 ( 26 , Lem. 7.2). Introduce 

(5.26) Bip, P ) = P M ^^^ Bip, P ). 

Then B(.,p) is a well defined meromorphic function, with divisor 

(5.27) iBi.,p)) = -Vp + Vfr-V p ,+V p , 
where v is defined via 

(5-28) a Eo iVo) = a Eo iV t ) + 2A Eo (p). 

Furthermore, 

(5.29) B( Poo ,p) = l, 

if 

6{z(poc,x,t) + 2A Eo (p)) 



(5.30) C p (x,t) 



Qizipvo,x,t)) 



Now, we can show how to conjugate the jump corresponding to one eigenvalue 
following |26j . 

Lemma 5.8. Assume that the Riemann-Hilbert problem for m has jump conditions 
near p and p* given by 

( 1 o\ 

m + (p) = rn_(p) 7 ( P ) 1 , p£ E s (p), 

(5.3i) y (p) - p \ 

I I 7(P*) \ 



m+ip) = m-ip) ^ Q <p)-p J , P G E e (p*), 
and satisfies a divisor condition 

(5.32) (mi)>-2V, (m 2 )>-^ A . 
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Then this Riemann-Hilbert problem is equivalent to a Riemann-Hilbert problem for 
fh which has jump conditions near p and p* given by 

' x B(p,p*){ix{p)-p) \ 
7(p)S(p*,p* 

,0 



m+ (p) = fh- (p) 



(5.33) 



m+(p) = m_(p) 



1 J 
1 0^ 

B(p" ,P*)(n(p)-p) 

l(P")B(p,p') J 



P G ^e(p), 
p£Z £ (p* 



divisor condition 

(5.34) (n*i)>-2V, (m 2 )>-% 

where T>c is defined via 

(5-35) a Bo p£) - " Bo (^ A) + 24 Bo (p). 

and all remaining data conjugated ( as in Lemma 1 5. 6\) by 

'B(p*,p*) 



(5.36) 



D(p) 



B(p,p* 

Proof. Denote by U the interior of Ti s (p). To turn 7 into 7 , introduce _D by 

'B(p*,p*) 



D(p) 



1 

7(p) 
jr(p)-p 





tt(p)-P n 
7(p) 



7<P* 



7r(p)-p 
1 



\ 7(P*) 

B(p*,p*) 
. B(p,p*) 



, Sfep* 


v 



else, 



and note that D(p) is meromorphic away from the two circles. Now set rh(p) = 
m(p)D(p). The claim about the divisors follows from noting, where the poles of 
B(p 7 p) are. □ 



Note that Lemma 15.81 can be applied iteratively to conjugate the eigenvalues 
Pj < C,{x/t): start with the poles p = p? and apply the lemma setting p — p\. This 
results in new poles p 1 = v. Then repeat this with p = p 1 , p = P2, and so on. 

All in all we will now make the following conjugation step: abbreviate 

, x, t) 



and introduce 



(5.37) D[p) = 



where 



lkip,x,t) 



1 



ULliPk ~ Pi) ^q{P*,X,t) 



Ik 



^(p,x,t) ^f^jDoip), l ^ k + l<e , Pk <C(x/t), 

■n\p)-Pk ! 

7 fc (p ,x,t) \ 

^~ Pk D (p), l<E , Pk < C(z/t), 







\lk (p*,x,t) 

[Doip), 

Doip) 



1 



else, 



T(p*,x,t) 
T(p,x,t) 
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Note that D(p) is meromorphic in K g \C(x/t) and that we have 
(5.38) D(p*)=(°, j^W/ 1 ' ' 



1 OJ ^' \l 0, 
Now we conjugate our problem using D(p): 

Theorem 5.9 (Conjugation). The function m 2 (p) — m(p)D(p), where D(p) is 
defined in (|5.37[) . is meromorphic away from C(x/t) and satisfies: 

(i) The jump condition 

(5.39) m\{p)=m 2 _{p)J 2 {p) 1 P eE, 
where the jump matrix is given by 

(5.40) J 2 (p) = D ^(p)- 1 J(p)D 0+ (p), 

(ii) the divisor conditions 

(5.41) (ro?) > -%*,*)., (n*a) > 

jumps corresponding to poles, except for possibly one if p k = C( x /t)> 
are exponentially decreasing. In that case we will keep the pole condition 
which is now of the form: 

lk{p,x,t) T(p*,x,t) 



m l P + — — — -rm 2 (p) \ >-V 0{x t y, near p k , 

/ 542 x V n(p)~p k T(p,x,t) ) 

f lk (p*,x,t) T(p,x,t) 2 2 \ 

(— n 7T m iW+ m 2(P)J >-Vt>( x , t ), near p k . 

\ TTlp) — Ok T(p*,X,t) ) - V 



"y k (p*,x,t) T{p,x,t) 
r(p)~ Pk T(p*,x,t)' 
(hi) i/ie symmetry condition 



m 2 (p*) = m 2 (p) 
(iv) and i/ie normalization 



1 

1 



m 2 ( Poo ) = (l 1). 

Proof. Invoking Lemma [5751 and f|4. 1 5[) we see that the jump matrix J 2 (p) is indeed 
given by (|5.40[) . The divisor conditions follow from the one for T(p, x, t) and m(p). 
Moreover, using Lemma 15.81 one easily sees that the jump corresponding to p k < 
£(x/t) (if any) is given by 

/ 1 T(p,x,t)t«(p)-Pk) \ 

J 2 ( P ) = I J 7*(P,»,t)r( P ^, t ) 1 , pe E e (p fc ), 
(5.43) ^ i " 7 

J 2 (p) = T(p%x,t)(7r(p)-p fc ) 1 , J3eS e (pfe), 

\ 7fc(p»,x,t)T(p,x,t) 

and by Lemma 15.61 the jump corresponding to p k > C,(x/t) (if any) reads 



(5.44) 



( 1 °\ 

■J (p) = afeWfj) 1 > P e s £ (p fc ), 

\T(p,x,t)(7r(p)-p fc ) i / 

/ 1 7 fc(p'V,*) , r(p^*) \ 

J 2 (p) = J T^, x ,t)U P )- Pk ) 1 p g S£(p *). 
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That is, all jumps corresponding to the poles pk ^ C(x/t) are exponentially decreas- 
ing. That the pole conditions arc of the form (|5.42p in the case pk — C{x/f) can 
be checked directly: just use the pole conditions of the original Riemann-Hilbert 
problem (|4. 1T[) and the divisor condition (|5.20[) for T(p,x,t). Furthermore, by 
(|4.18l) and (|5.38[) one checks that the symmetry condition for m 2 is fulfilled. From 
T(poo, x,t) = 1 we finally deduce 

(5.45) m 2 { Poo ) = m(poo) = (l l) , 

which finishes the proof. □ 
Forp e E\C(x/*) = EnTr-^^x/^.oo)) the jump matrix J can be factorized 

as 

J 2 = (b-)- l b+, 

where b± = D^ 1 b±DQ, that is, 
(5.46) 

- (l gteL iz(p*)e(p»)e-«»W\ - / 1 o\ 

L > I J- 4+ ={m»^)^ 'J- 

For p € C(ir/f) = E n tt" 1 ((-oo, C(^A))) 

we can factorize J 2 in the following way 

J 2 = (B-)~ 1 B + , 

where B± = D^. 1 B±D±, that is, 
(5.47) 

/ 1 (A ~ /l T + (p,x,t) R(p*)<5>(P*) r -t<t>(v)\ 

B. = [ T-y,,,t)fl(p)e(p) c «^ J, B + = J t +M e 



Note that by T(p, x, t) = T(p, x, t) we have 



KifC\ T-(p*,x,t) _ T-(p*,x,t) 1 T+(p,x,t) 

[ ' T + (p,x,t) T-(p,x,t) 1-\R( P )\ 2 T + (p,x,t)> P eU ^' 

respectively 



(5.49) T+jP'*'*) T + fp,*,t) 1 T-(p*,z,t) c(x/t) _ 

1 ' T.(p*,x,t) T + (p*,x,t)l~\R(p)\ 2 T_(p*,x,tY P Kl) 

We are now able to redefine the Riemann-Hilbert problem for m (p) in such a 
way that the jumps of the new Riemann-Hilbert problem will lie in the regions 
where they are exponentially close to the identity for large times. The following 
theorem can be checked by straightforward calculations: 

Theorem 5.10 (Deformation). Define m 3 (p) by 



(5.50) m 3 (p) = < 



m 2 {p)B + {p)-\ pEDkU D ju k < j, 

m 2 {p)B_{p)-\ P eD* k U D* x , k < j, 

m 2 (p)b+{p)- 1 , peD k U D j2 , k > j, 

)* 

f2 



m 2 {p)b_{p)-\ peD* k UD* j2 , k>j, 
m 2 (p), else, 



where the matrices b± and B± are defined in (|5.46j) and (|5.47[) . respectively. Here 
we assume that the deformed contour is sufficiently close to the original one. Then 
the function m 3 (p) satisfies: 
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(i) The jump condition 

(5.51) m+(p) = m 3 _(p)J 3 (p), for p <G S, 

where the jump matrix J 3 is given by 

f B+(p), P e C k UC jU k < j, 
B-(p)-\ peCluq,, k<j, 
b+(p), peC k UC j2 , k> j, 

6_(p)-\ P eC£uc; 2 , k>j, 

J 2 (p), else, 

(ii) the divisor conditions 

(5.53) (m?) > (m 3 ,) > -% (x , t) , 

J7ie jumps on the small circles around the eigenvalues remain unchanged. 

(iii) The symmetry condition 



(5.52) J 3 {p) = < 



1 

1 



(5.54) m 3 (p*) = m 3 (p) 

(iv) and the normalization 

(5.55) m 3 ( Poo ) = (1 1) . 

Here we have assumed that the reflection coefficient R(p) appearing in the jump 
matrices admits an analytic extension to the corresponding regions. Of course this 
is not true in general, but we can always evade this obstacle by approximating R{p) 
by analytic functions. We relegate the details to Section [7] 

The crucial observation now is that the jumps J 3 on the oriented paths C k , C£ 
are of the form I + exponentially small asymptotically as t — > oo, at least away 
from the stationary phase points Zj, z* . We thus hope we can simply replace these 
jumps by the identity matrix (asymptotically as t — ¥ oo) implying that the solution 
should asymptotically be given by the constant vector (l l) . That this can in fact 
be done will be shown in the next section by explicitly computing the contribution 
of the stationary phase points thereby showing that they are of the order O^ 1 / 2 ), 
that is, 

m 3 (p) = (l l) + 0{t- 1 ' 2 ) 

uniformly for p away from the jump contour. Hence all which remains to obtain 
the leading term Vi in Theorem 12.21 is to trace back the definitions of m 3 and m 2 
and comparing with (|4. 13[) . First of all, since m 3 and m 2 coincide near p^ we have 

m 2 ( P ) = (i i)+o(r 1 / 2 ) 

uniformly for p in a neighborhood of Poo. Consequently, by the definition of m 2 
(see Theorem 15. 9p . we have 

m(p) = (Tip^x^t)- 1 Tip.x^-^+Oir 1 ' 2 ) 

again uniformly for p in a neighborhood of poo- Finally, using the expansion of 
T(p, x, t) near p^ (see Lemma [S~5|) and then comparing the last identity with (|4.13[) 
shows 

/>oo 

(5.56) / (V-V q )(y,t)dy = 2iT 1 (x,t)+0(t- 1 ^ 2 ), 
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where T\ is defined via (|5.22|) . that is, 

T 1 (x,t)=- 2 T ^,0 + ^- f log(l-|i?| 2 V Poo , 

Pk <ax/t) Je ^) 2m J c(x/t) 



id T In 



0(z(Poo,x,t)) 



Similarly one obtains 

(5.57) (V-V q )(x,t) = 0(t- 1 / 2 ), 

by using P~2"7]) instead of P~T3|) . 

Hence we have proven the leading term in Theorem 12.21 the next term will be 
computed in Section [B] 

Case (ii): The soliton region. In the case where no stationary phase points lie 
in the spectrum the situation is similar to the case (i). In fact, it is much simpler 
since there is no contribution from the stationary phase points: There is a gap (the 
j-th gap, say) in which two stationary phase points exist. Similarly as in case (i) 
an investigation of the sign of Re((/>) shows the following: 

fRe(0(p)) > 0, P eD k , k <j, 
\Re(<b( P )) < 0, P eD k ,k> j. 

Now we construct "lens-type" contours C k (as shown in Figure [2|) around every 
single band lying to the left of the j-th gap and make use of the factorization 
J 2 = (6_) _1 6+, where the matrices 6_ and 6+ are defined in (|5.46[) . We also 
construct such "lens-type" contours C k around every single band lying to the right 
of the j-th gap and make use of the factorization J 2 = (B-)~ 1 B+ with the matrices 
£?_ and B + given by Q5.47p . Indeed, in place of (I5.50|) we set 

m 2 (p)B- 1 {p), peD k , k< j, 
m 2 (p)BZ 1 {p), peD* k , k<j, 



(5.58) m 3 (p) 



'ki 

m 2 (p)b+ 1 (p) 1 p£ D k , k> j, 
m 2 {p)bZ 1 {p), p£ D* k , k> j, 

,2/ 



m (p), else. 

Now we are ready to prove Thcorcm l2.Il bv applying Thcorcm lA.8l in the following 
way: 

If \C,{x/t) — pk\ > £ for all k we can choose 7q = and Wq by removing all 
jumps corresponding to poles from w . The error between the solutions of and 
Wq is exponentially small in the sense of Theorem IA.8) that is, — WqIIoo < 
0(t ) for any I > 1. We have the one soliton solution (cf. Lemma f4.6|) rho(p) = 
(f(p*,x,t) f(p,x,t)), where f(p) = 1 for p large enough. Using Lemma [5.51 we 
compute 

[Tip^x,^- 1 
V T(p,x,ty 



m(p) =m (p) t n Ti ^ , M j 



= ( l + IlMl + 0(z- 1 ) l~ll^l+0(z- 
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Comparing this expression with (|4.13p yields 

(V ~ V q )(y, t)dy = 2iT 1 (x, t) + 0(r n ), 
and thus by our definition of the limiting solution we finally have 

V - Vi)(y,t)dy = [°° f(V-V q )(y,t) - (V t - V q )(y,t))dy = 0(t""), 

J X 

for any n > 1 if R(p) has an analytic extension. This proves the second part of the 
theorem. 

If \((x/t) — pk\ < e for some k, we choose 7q = % and Wq = 0. Again we 
conclude that the error between the solutions of w l and Wq is exponentially small, 
that is, — WqWoo < 0(t~ l ), for any I > 1. By Lemma FO)l we have the one soliton 
solution m (p) = (f(p*,x,t) f(p,x,t)), with 

i, ^ - , Ik ^Lc k (Pk,x,t)W {Xit) (^ Ck (p k ,x,t),iP ltCk (p,x,t)) 

f(P,X,t) = 1 + : — , 

z-pk ipi,c k {p,x,t)ci,k{x,t) 
for p large enough, where jk is defined as in (|2.13[) . We will again use 

™(rt^)( r(p *r r ' T(Pi • r .) = tea). 

and now expand f(p) as in the proof of Lemma 14.61 Finally a comparison with 
(|4~T3|) yields 

V - V q )(y, t)dy = 2m( X , t) + 2 ^ l ^ k ^ t)2 + 0{t- n ), 

ci,k{x,t) 

and hence by our definition of the limiting solution (|2.8[) we obtain (I2.10P for any 
n > 1 if R(p) has an analytic extension. Similarly we obtain (I2.11j) by using (|4.27[) 
instead of (|4.13p . 

Case (iii): The transitional region. In the special case where the two stationary 
phase points coincide (so Zj — z* — Ek for some fc) the Riemann-Hilbert problem 
arising above is of a different nature. We expect a similar behavior as in the constant 
background case [10] (cf. also the discussion in [22] )■ However, we will not treat 
this case here. 

6. The "local" Riemann-Hilbert problems on the small crosses 

In the previous section we have reduced everything to the solution of the Riemann- 
Hilbert problem 

m 3 + (p) = m 3 _(p)J 3 (p), 

,3\ 



K)>-%x,t)«, (m 2 ) > -T>^ x ,t), 

1 N 

1 o, 



m (j>*) = m (p) 
m 3 ( Poo ) = (l 1), 



where the jump matrix J 3 is given by (|5 . 52[) . We have performed a deformation in 
such a way that the jumps J 3 on the oriented paths Ct, for k ^ j are of the 
form "I + exponentially small" asymptotically as t — !• oo. The same is true for the 
oriented paths Cj%, Cj2, C* 1: C* 2 at least away from the stationary phase points 



30 



A. MIKIKITS-LEITNER AND G. TESCHL 





Figure 3. The small cross containing the stationary phase point 
Zj and its flipping image containing z*. Views from the top and 
bottom sheet. Dotted curves lie in the bottom sheet. 

The purpose of this section will be to derive the actual asymptotic rate at 
which m 3 (p) — > (l l) following again [22]. The jump contour near the stationary 
phase points (cf. Figured]) will be denoted by Y, c (zj) and On these crosses 

the jumps read 

(\ -Z_ iTO* -t<t>\ 



J 3 

3.1 



J 3 — -S- 1 — ( T* R& t cf> i j ) P e 7 



1 
-i?ee*^ 1 



o i ^ c ^-2- 

To reduce our Riemann-Hilbert problem to the one corresponding to the two 
crosses we proceed as follows: We take a small disc D around Zj{x/t) and project it 
to the complex plane using the canonical projection 7r. Now consider the (holomor- 
phic) Riemann-Hilbert problem in the complex plane with the very jump obtained 
by projection and normalize it to be I near oo. 

The corresponding Riemann-Hilbert problem is solved in [251 Appendix A]. To 
apply [25j Theorem A.l] we need the behavior of the jump matrix J 3 , that is, the 
behavior of T(p,x,t) near the stationary phase points Zj and z*. 

The following lemma gives more information on the singularities of T(p, x, t) near 
the stationary phase points Zj, j = 0, . . . , g and the band edges Ej, j = 0, . . . , 2g+ 1 
(setting E 2g +i = oo). 

Lemma 6.1. For p near a stationary phase point Zj or z* (not equal to a band 
edge) we have 

(6.2) T{p,x,t) = {z-z.i) ± ' w e ± {z) 1 p=(z,±), 
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where e (z) has continuous limits near Zj and 

(6.3) ^=-^log(l-|i?( 2j )| 2 )>0. 

Z7T 

Here (z — Zj) " = exp(±i^log(z — Zj)), where the branch cut of the logarithm is 
along the negative real axis. 

For p near a band edge Ek € C(xjt) we have 

(6.4) T(p,x,t)=T ± \z)e ± (z), p=(z,±), 

where (z) is holomorphic near Ek if none of the Vj is equal to Ek and e± (z) has 
a first order pole at Ek = Vj else. 

Proof. By factorizing the jump according to 1 — |i?(p)| 2 = (1 — | R ( z j ) | 2 ) ^ ~| 1 2 
one can rewrite (|5.19p . Then one considers the Abelian differential w pp » for p € 
lC g \{poc} which is explicitly given by the formula (I3.21j) . One has 

(6.5) \ f uj pp , = ±log(z - zj) ± a(zj) + 0(z - Zj), p=(z,±) 

1 JC(x/t) 

and thus 

(6.6) / u PPoo = ±log(z- Zj)±a(zj) + 0(z- zj), p=(z,±), 

JC(x/t) 

since / c(x/t) fu PPoo = \ J c{x/t) fuj pp * for any symmetric function f(q) = f(q*). 
From this the first claim follows. For the second claim note that the function 



tip) 



T(z), p = (z,+)en+, 
T{z)-\ p=(«,-)en_, 

satisfies the following (holomorphic) Riemann-Hilbert problem 

t + {p)=t^{p)(i-\R{ P )\ 2 ) 1 peE, 

t{Poo) = 1- 

Thus T(p,x,t)/t(p) has no jump along C(x/t) and is therefore holomorphic near 
C(x/t) away from band edges Ek = Vj (where there is a simple pole) by the Schwarz 
reflection principle. □ 

Moreover, 

Lemma 6.2. We have 

(6.7) e ± (z) =W{z), p = (z,±)€'E\C(x/t) 

and 
(6.8) 

z(poo,x,t) + S(x/t)) 9(z(Zj,X,t)) 



-(zj) = cxp (wa(zj))- 



U(Poo,a;,i)) 0(z(zj,x,t) + 8(x/t)) 



• exp ( - £ I1L u « *> + ^ L /t) log (t^Sf)^ 



p k <C(x/t) 



where 



(6.9) a(Zj) = lim (i [ w pp . - log (tt(p) - Zj)). 

P^Zj J C (x/t) ' 

Here ct{zj) £ R and w pp * is real on C{x/t). 
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Proof. The first claim follows from the fact that 



T(p*,x,t) =T(jp,x,t) =T(p,x,t) forpe Y,\C(x/t). 

For the second claim just follow the argument used in the proof of the previous 
lemma. □ 

By Lemma |6. II one deduces that near the stationary phase points the jumps are 
given by 

1 0\ 



(6.10) 



1 ON 

-2ii/ 



^l{ z - Zj) y re** l) ' Peij2 ' 



,0 



2ir/ 

J2' 



where (cf. eq. (16.2 



(6.11) r = R( Z] M Zj , x, tf^l " 

The error terms will satisfy appropriate Holder estimates, that is 

(6.12) ||B+(p)-B + (p)|| <C|z-z,r, p=(z,+)eC jl , 

for any a < 1 and similarly for the other matrices. Thus the assumptions of [231 
Theorem A.l] are satisfied and we can conclude that the solution on 7r(S c (z :) )) is 
of the form 

(6.13) M(z)=I+AJ_ + o(r a ), 



: - zj tV2 



where 



m o 



(6.14) Mo = 

(6.15) 0(t) - v ^i(T/4-«s(r)+«r«(r(U/))) e -t*(^) t -i Vj 

and 1/2 < a < 1. Now we lift this solution in the complex plane back to the small 
disc D on the Riemann surface JC g by setting 

f M(z), peD, 

( " 6) M(p) = \ (! i) *<«>(! J). 

Thus we conclude that the solution on E (zj) is given by 
(6.17) tfC W = I V^ + °( ra '' P =(*,+). 
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(6.19) m 4 (p) = < 



and the one on S c (z*) reads 

(6.18) MC {p)=l+ _LJ^_ + 0{t - a): p =(*,_). 

Then 

' m 3 (p)M c (p)~ 1 , peD, 
m 3 (p)M c {p)- 1 , peD*, 
m 3 (p), else, 

has no jump inside DUD* but jumps on the boundary 3D U dD*. All jumps 
outside D U D* are of the form I + exponentially small and the jump on dD U dD* 
is of the form I + 0(i -1 / 2 ). In order to identify the leading behavior it remains 
to rewrite the Riemann-Hilbert problem for m 4 as a singular integral equation 
following SectionEl Let the operator C w i : L 2 (£ 4 ) — > L 2 (E 4 ) be defined by 

(6.20) C w ,f = C^(fw i ) 

for a vector valued /, where u> 4 = J 4 — I and 

if - fn- ,Poa o \ 



3.21) (C±/)(g)= Um / /Of, 0|: 



fi**» 



are the Cauchy operators for our Riemann surface. In particular, Qp' q is the Cauchy 
kernel given by 

(6-22) ^=^ + ]T/f>)G, 

i=i 

where 

(6.23) #<(p) = £<*(*) r^ i0 . 

Here w 9i o is the (normalized) Abelian differential of the second kind with a second 
order pole at q (cf. Remark 16.41 below) and ui pq denotes the Abelian differential of 
the third kind with simple poles at p and q. Note that Ij' 9 (p) has first order poles 
at the points P. 

The constants Cji(v) are chosen such that f2fj' 9 is single valued, that is, 

-l 



dm 



(6.24) £ c *0')-^7a , 

where Cfe(j) are defined in (|3.6|) (cf. Lemma |A~3|) . 

Consider the solution fx of the singular integral equation 

(6.25) M=(l l)+C w 4/x in L 2 (£ 4 ). 
Then the solution of our Riemann-Hilbert problem is given by 

(6.26) m\p) = (l l^-L^yu, 4 ^. 

By ||w 4 ||oo = 0(tT x l 2 ) Neumann's formula implies 

(6.27) f x\q) = (I-C w ,)- 1 (l 1) = (1 l)+0{t- 1 ' 2 ). 
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Moreover, 

f_jv|o — i_ + (t- a ), pedD 



Z-Zj t 1 / 2 

-^-^ + 0(<- Q ), pedD*. 



(6.28) w 4 (p) = 

Hence we obtain 

t 1 / 2 2tti J 9d K-Zj p 

(1 l)Mn 1 /" 1 



tVa 2niJ dD ,n-z J -v + U[t ' 



= (i i)-^#0l(^-£#^) + o(*--) 

x (i0fi£*' p ~ - i/Jftf^ (z,-) + i£fif' Po ° (z*)) 

(6.29) + 0(t~ a ). 

Since we need the asymptotic expansions around p^ we note 
Lemma 6.3. We have 

(6.30) O^-te) - Af(z,)C + Af(z,)C 2 + 0(C 3 ) 

/or C = z -1 / 2 femj £/ie Zocc^ chart near p^ and 

a a 

(6.31) Ay(z,-) = w PDO! o(^) - ^^c M (^)a g _i(^)Cfe(zj), 

1 5 9 

(6.32) Af(z 3 ) = Wp^ifo) - - ^^cfe^^)^^), 

k=l 1=1 

where 0J q .k, k = 0, 1, . . . , is an Abelian differential of the second kind with a single 
pole of order k + 2 ai q and a g -i{i>t) denotes a constant defined in Remark \6.f\ 
below. 

Proof. Use the local coordinate Q = z -1 / 2 near p^ — (oo,oo) and expand the 
differential u> PPoo like it is done in |35( Theorem 4.1] and J wp fj o by using the 
expression (|6.33[) . For g > 1 one gets 

^,o(C) = -« 9 -i(^)-C + o(C 2 ) 

and thus the claimed formulas for Aj(zj) and Aj(zj) follow. □ 

Remark 6.4. The Abelian differential appearing in the previous lemma is explicitly 
given by 

<+x + RlilM + -") + p ^o ■ (tt - V? 

(6.33) Looa = — 775 dir, 
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with Ppfi a polynomial of degree g—1 which has to be determined from the normal- 
ization. We will use the notation 

9-1 

(6.34) P fij0 (z) = J2<*jP> j - 

3=0 

Concerning the Abelian differential w Poo ,o we refer to (|3.15[) . The differential uJ Poa ,i 
is given by 

(6-35) ^, 1= (_:^ + P W ) *L_, 

H 2g+1 

where P Poo ,i is a polynomial of degree g—1 which has to be determined by the 
vanishing a j -periods as usual. 



Note that the following relations are valid 
(6.36) 



and 

(6.37) c k e(y*) = -Ckt(0, Ck(zj) = ~Ck(zj)- 

Moreover, the coefficients a 3 - (/>), j = 0, ...,<? — 1 of the polynomial Po$ fulfill the 
relation 

(6.38) aj (P*) = -ajp), j g 1. 

Now we come to prove Theorem l2.2l As in the previous section, the asymptotics 
can be read off by using 



l 







(6.39) ro(p) = m\p) T ^ 

for p near poo and comparing with (|4.13|) . From that one deduces 
m 2 (p) =TO2(p) t (p) _1 = 



where we have used (I6.29[) . (|6.30[) and (j5.21j) . Comparing this asymptotic expansion 
with (|4~T3f yields 



+ °V - V q )(y)dy = 2 y ^L_(/3Affe) - (3Af(z*)) + 2\T x {x,t) + 0(^ Q ). 
Invoking (j(H6l) . (f6T37|) and (f6T38|) one gets 
(,40, 4 (z -) = -Af fe ), 

Af {Zj) = At(z,), 

and therefore 

1 



3.41 



{V-V q ){y,t)dy=U-^— t Re((3(x,t))Af(x > t) + 2iT 1 (x,t) + 0(t- a ). 



3(3 
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Finally, using the definition of the limiting solution (|2.8j) proves the claim. Note 
that one obtains the same result if one compares the expressions for the component 
mi. 

Similarly one obtains (|2.17|) by using (|4.27l) instead of (|4.13j) . 

7. Analytic Approximation 

In this section we want to show how to get rid of the analyticity assumption 
on the reflection coefficient R(p). To this end we will split R{p) into an analytic 
part J? ajt plus a small rest R r t following the ideas of [5] (see also Sect. 6]). 
The analytic part will be moved to regions of the Riemann surface while the rest 
remains on £ = 7r _1 (cr(i/ 9 )) . This needs to be done in such a way that the rest is 
of 0(i _1 ) and the growth of the analytic part can be controlled by the decay of the 
phase. 

In order to avoid problems when one of the poles Vj hits E, we have to make the 
approximation in such a way that the nonanalytic rest vanishes at the band edges. 
That is, split R according to 

m =fl(^) ir^ )( fHr~-? + {z " E2j+l){E2j ~ i} 



(E 2j+ i - E 2j )(z - i) v '™ (E 2j - E 2j+1 )(z - i) 

(7.1) +Q J (p)R( P ), p=(z,±), 

where Qj(p) is a rational function with hrst order zeros at E 2 j 7 E 2 j + i and with all 
other zeros and poles away from S, and approximate R. Note that if R G C'(£), 
then R G C'-^S). 

We will use different splittings for different bands depending on whether the 
band contains our stationary phase point Zj(x/t) or not. We will begin with some 
preparatory lemmas. 

For the bands containing no stationary phase points we will use a splitting based 
on the following Fourier transform associated with the background operator H q . 
Given R G C l (£) we can write 

(7.2) R(p) = [ R{x)^ q (p,x,0)dx, 

JR 

where ip q (p,x,t) denotes the time-dependent Baker- Akhiezer function and (cf. [14] . 

M) 

i r in 9 iMp)-a*7\) 

(7.3) R(x) = —& R(p)%(p,x,Q) LL > =1 \ W / 3; <fer(p). 

2tti h 2R 1 2 / g 2 +l (p) 

If we make use of p,14p . the above expression for R(p) is of the form 

(7.4) R(p)= / R(x)9 q (p,x,0)e^p(ixk(p))dx. 



where k(p) = — i J^ w Poc ,o an d O q (p 7 n,t) collects the remaining parts in p,14[) . 
Using k(p) as a new coordinate and performing I integration by parts one obtains 

(cf. m) 

,„ ,A, . const 

(7-5) \R(x)\ < 

provided R G C l (E). 



ASYMPTOTICS OF PERTURBED FINITE-GAP KDV SOLUTIONS 



37 



Lemma 7.1. Suppose R e ^(M), x l R{x) G L 1 (R) and let /3 > be given. Then 
we can split R(p) according to 

R(p)=Ra,t(p)+Rr,t(p), 

such that R a ,t(p) is analytic for in the region < Im(fc(p)) < e and 

(7.6) |i? a , 4 (p)e-' 3t | = 0(t- 1 ), < lm(k(p)) < e, 

(7.7) \R r . t {p)\=0(tr l ), peS. 
Proof. We choose 

/>OG 

R a ,t(p)= / R{x)O q {p, x, 0)exp(ixk(p))dx 

Jx=~K{t) 

with if(t) = for some positive /?o < /?• Then, for < Im(fc(p)) < e, 

/•OO 

|i?a,t(fc)e-' 3t | < Ce"^ / |i?(x)|e- Im(fc(p))x dx 

<Ce-P t e K ( t >\\F\\i=\\R\\ 1 e-V-M t , 
which proves the first claim. Similarly, for p 6 S, 

ip (v\\<rr ^L^ll w ^^ ll^(-^)lli ^ g 

□ 

For the band which contains Zj{x/t) we need to take the small vicinities of the 
stationary phase points into account. Since the phase is cubic near these points, we 
cannot use it to dominate the exponential growth of the analytic part away from S. 
Hence we will take the phase as a new variable and use the Fourier transform with 
respect to this new variable. Since this change of coordinates is singular near the 
stationary phase points, there is a price we have to pay, namely, requiring additional 
smoothness for R{p). 

Without loss of generality we will choose the path of integration in our phase 
<fr(p), defined in (I4.21[) , such that 4>(p) is continuous (and thus analytic) in Dj t i 
with continuous limits on the boundary (cf. Figure [T]). We begin with 

Lemma 7.2. Suppose R(p) G C J (S). Then we can split R(p) according to 

(7.8) R(p) = Ro(p) + ^i-H(p), p=( 2 ,±)eSnD j , 1 , 

where Ro(p) is a real rational function on M such that H(p) vanishes at Zj, z* of 
order three and has a Fourier transform 

(7.9) H(p) = [ H(x)e xcl>{p) dx, 

Js. 

with xH(x) integrable. Here <f> denotes the phase defined in (|4.21[) . 

1/2 

Proof. We begin by choosing a rational function Ro{p) = a{z) + b(z)R 2g+1 (p) with 
p = (z, ±) such that a(z), b{z) are real- valued rational functions which are chosen 
such that a(z) matches the values of Re(i?(p)) and its first four derivatives at Zj 
and ]r 1 b{z)R\^ +1 {p) matches the values of Jm(R(p)) and its first four derivatives 
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at Zj. Moreover, all poles are chosen away from S. Since R(p) is C 5 we infer that 
H(p) £ C 4 (E) and it vanishes together with its first three derivatives at Zj, z* . 

Note that <f)(p)/i, where <fr is defined in (|4.21[) has a maximum at z* and a mini- 
mum at Zj. Thus the phase <f>(p)/i restricted to EH-D^i gives a one to one coordinate 
transform £ (1 Dj t i — > [<j)(z*)/i, <f>(zj)/i] and we can hence express H(p) in this new 
coordinate (setting it equal to zero outside this interval) . The coordinate transform 
locally looks like a cube root near Zj and z*, however, due to our assumption that 
H vanishes there, H is still C 2 in this new coordinate and the Fourier transform 
with respect to this new coordinates exists and has the required properties. □ 

Moreover, as in Lemma 17. II we obtain: 

Lemma 7.3. Let H(p) be as in the previous lemma. Then we can split H(p) 
according to H(p) = H a ^(p) + H r ^(p) such that H a ^(p) is analytic in the region 
Re(0(p)) < and 

(7.10) \H a , t (p)e^ 2 \ = 0(1), peTJJJ, \H r , t (p)\ - Oir 1 ), P eE. 

Proof. We choose H a>t (p) = $™ = _ K{t) H (x^^ dx with K{t) = t/2. Then we can 
proceed as in Lemma 17.11 

|#a,tb)e 0(p)t/2 | < ||.ff|| 1 |e-- K W*W+*(i»)t/2| < 

and 

1 f°° " C 

\H r M\ < wTT / X\H(-X)\dx< -. 

K{t) J x= K(t) t 

□ 



Clearly an analogous splitting exists for p e En D j2 • 

Now we are ready for our analytic approximation step. First of all recall that our 
jump is given in terms b± and B± defined in (|5.14[) and (|5.16p . respectively. While 
b± are already in the correct form for our purpose, this is not true for B± since 
they contain the non-analytic expression |T(p)| 2 . To remedy this we will rewrite 
B± in terms of the left rather than the right scattering data. For this purpose 
let us use the notation R r (p) = R+(p) for the right and Ri{p) = R-(p) for the 
left reflection coefficient. Moreover, let T r (p,x,t) = T(p,x,t) be the right and 
Ti(p,x,t) = T(p)/T r (jp,x,t) be the left partial transmission coefficient. 

With this notation we have 



(7.11) J 2 ( P ) = 

where 



b-^b+ip), 7T(p) > ((x/t), 

B-(p)- l B + (p), n(p)<C(x/t), 
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and 



B_ = 
B+ = 

Using (|4.8p we can write 

B_ = 

B, = 



T r ,-(j>',x,t) R r (p)e(p) „t 0(p) 
T ri .(p,x,t) \T(p)\* e 

2 e 



T r ,+ (p* .a;,*) |T(p)| 
1 



*<Mp) 



v Ti{p.x,t) 

1 



Now we split i? r (p) = R a ,t{p)+Rr,t{p) by splitting -R r (p) defined via (|7.1I) according 
to Lemma 1 7. II for 7r(p) G [i?2fc, -Esfc+i] with k < j (i.e., not containing ((x/t)) 
and according to Lemma [7.31 for tt(p) £ \E2j,C,{x/t)\. In the same way we split 
Ri{p) = Ra,t{p) + Rr,t{p) for i"(p) G [C(i/t),^j+i] and 7r(p) G [£ 2fe ,£ 2 fc + i] with 
k > j. For /3 in Lemma 17.11 we can choose 



(7.12) 



_ I mm peCfc 



-Re(0(p)) > 0, tt(p) > C(a:/t), 

7T(p) < C( X /t). 



min p6Cfc Re((j)(p)) > 0, 
In this way we obtain 

b±(p) = b a ,t,±{p)b r ,t,±(p) = br,t,±(p)b a ,t,±(p), 
B±( P ) = B a .t,±{p) B r,t.±{p) = Br <t ,±(P)Ba,t,±(P)- 

Here b a ,t,±{p), b r ,t,±(p) denote the matrices obtained from b±(p) by replacing R r {p) 
with Rajip), R r ,t(p), respectively. Similarly for -B a ,t,±(p); Br,t,±(p)- Now we can 
move the analytic parts into regions of the Riemann surface as in Section [5] while 
leaving the rest on E. Hence, rather than (|5.52j) . the jump now reads 



(7.13) 



J 3 (P) = { 



ba,t,+ (p), 

ba,t,-{P)~ l 
6r,t,-(p)~ 1 & 

B a ,t,+ (p), 

Ba^-ip)- 1 



r,t,+ 



(P), 



P G C fc , tt(p) > C(a;/«), 
peQ, tt(p) > C(i/t), 
PGtt-^^x/^.+oo)), 
p G C fe , ?r(p) < C(x/t), 
P eC* k , nip) < C(x/t), 



S r , t ,_(p)- 1 J B r , t , + (p), pGTT-H-oo.CCx/*))). 



By construction R a ,t{p) = Ro{p) + (tt(p) — ^{zj))H a ,t (p) will satisfy the required 
Lipschitz estimate in a vicinity of the stationary phase points (uniformly in t) and 
the jump will be J 3 (p) = I + 0(t _1 ). The remaining parts of £ can be handled 
analogously and hence we can proceed as in Section [51 



Appendix A. Singular integral equations 

In the complex plane, the solution of a Riemann-Hilbert problem can be reduced 
to the solution of a singular integral equation (see [3]). In our case the underlying 
space is a Riemann surface M. The purpose of this appendix is to generalize this 
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approach to meromorphic vector Riemann-Hilbert problem with simple poles at p, 
p* of the type 

m + (p) — m-(p)J(p), p G E, 

(mi) > —T>n* - V p , (m 2 ) > -Vp, - V p , , 



) 

,1/2 



nLi (P - Vk)n(p)- P % (P* ) ' 
'0 1 N 
.1 o, 



mi(p) + m 2 (p) > -Vp,, near p*, 



771 (p*) = 771 (p) 



m(poo) = (l l) , 

Since we require the symmetry condition (|4.18p for our Riemann-Hilbert problems, 
we need to adapt the usual Cauchy kernel to preserve this symmetry. Moreover, 
we keep the single soliton as an inhomogeneous term which will play the role of the 
leading asymptotics in our applications. 

Concerning the jump contour E and the jump matrix J we will make the following 
assumptions: 

Hypothesis H. A.l. Let E consist of a finite number of smooth oriented finite 
curves in M which intersect at most finitely many times with all intersections being 
transversal. Assume that the contour E does not contain any of the points fi and 
is invariant under p i-> p* . It is oriented such that under the mapping p H> p* se- 
quences converging from the positive sided to E are mapped to sequences converging 
to the negative side. The divisor Vp is nonspecial. 

The jump matrix J is nonsingular and can be factorized according to J = 
bZ_ l b + = (I — ?«_) _1 (I + w + ), where w± = ±(6± — 1) are Holder continuous and 
satisfy 



(A.2) ,r U n (J JW(p)(5 J), P^. 



Moreover, 

(A.3) IMloo = ll^+IU°°(£) + ||tt>-||£°°(E) < 0°- 

Remark A.2. The assumption that that none of the poles p lie on our contour E 
can be made without loss of generality if the jump is analytic since we can move the 
contour a little without changing the value at poo (which is the only value we are 
eventually interested in). Alternatively, the case where one (or more) of the poles 
fjj lies on E can be included if one assumes that w± has a first order zero at pbj . 
In fact, in this case one can replace p,(s) by p(s) = (tt(s) — pj)p(s) and w±(s) by 
w±(s) = (ir(s) - p J )~ 1 w ± (s). 

Otherwise one could also assume that the matrices w± are Holder continuous and 
vanish at such points. Then one can work with the weighted measure — ii?2g+i(p)d7r 
on E. In fact, one can show that the Cauchy operators are still bounded in this 
weighted Hilbert space (cf. |19[ Thm. 4.1],). 
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Our first step is to replace the classical Cauchy kernel by a "generalized" Cauchy 
kernel appropriate to our Riemann surface. In order to get a single valued kernel 
we need again to admit g poles. We follow the construction from (3UJ Sec. 4]/ 

Lemma A. 3 ( 22, 25 ). Let be nonspecial and introduce the differential 

(A.4) ^ P = Wpp + ]T/f P (p)G, 

where 



(A.5) /f>) = f 



Here w 9i o is the (normalized) Abelian differential of the second kind with a second 
order pole at q (cf. Remark \6.4\ l and the matrix Cjg is defined as the inverse matrix 
of r\i(Jij), where Q = r\n(z)dz is the chart expression in a local chart near jij (the 
same chart used to define wa.^j. 



Then Slp' P is single valued as a function of p with first order poles at the points 




Next we show that the Cauchy kernel introduced in (|A.4[) has indeed the correct 
properties. We will abbreviate U>(E) = LP(E,C 2 ). 

Theorem A.4 ([22"1 125]). Set 

(a.6) of' p = r 1 ^" 

and define the matrix operators as follows. Given a 2 x 2 matrix f defined on S 
with Holder continuous entries, let 

(A.7) (Cf)(p) = ±;Jj(s)ti' P , for p^E, 

and 

(A.8) (C±/)(?)= lim (Cf)(p) 

from the left and right ofT, respectively (with respect to its orientation). Then 

(i) The operators C± are given by the Plemelj formulas 

(C+f)(q) - = f(q), 

( C+ f)(q) + (C-f){q) = -l fnl' P , 

and extend to bounded operators on L 2 (E). Here £ denotes the principal 
value integral, as usual, and the bound can be chosen independent of the 
divisor as long as it stays some finite distance away from E. 

(ii) Cf is a meromorphic function off E, with divisor given by ((Cf)ji) > 
-2V and{{Cf) j2 )>-V t 

(iii) (Cf)( P *) = (0 *), (Cf){p) = (* 0). 
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Furthermore, C restricts to L 2 (T,), that is 

(A.9) (Cf)(p*) = (Cf)(p)(° JY P e 

for f £ L 2 (T,) and if w± satisfy (|A.2[) we also have 



1 




pes. 



(A.lO) C ± (fw T )( P *) = C T (fw ± )( P ) 

Now, let the operator C w : L 2 S (Y>) -> L 2 (T,) be defined by 
(A.ll) C w f = C+{fw-)+C-(fw+), /€L2(E), 

for a 2 x 2 matrix valued /, where 

w+ = o+ — I and w- = I — 6- . 
Recall from Lemma 14.61 that the unique solution corresponding to J = I is given by 

mofr) = (/(?*) /(?)), 
for some given f(p) with (/) > —"Dp, — T> p * . Since we assumed T>^ to be away from 
S, we clearly have mo G 

Theorem A. 5 ([22l [25]). Assume Hvvothesis \A.l\ and let mo € C 2 6e given. 
Suppose m solves the Riemann-Hilbert problem (|A.lj) . T/ien 

(A.12) m(p) = (1 - co)mo(p) + ^ f n(s)(w + (s) + W-{s))Q§ 



ii. p 



wh 



ere 



fi = m+b , 1 = m_6_ 1 and Co = (7^—7 / l^(s)(w+(s) + w-(s))Q^) . 

-Here (m) j denotes the j 'th component of a vector. Furthermore, il solves 

(A.13) (l-0> = (l-co)m (p). 

Conversely, suppose fx solves 
(A.14) (I-C,„)A = mo, 

and 

5 - jf A(*)(tfl+w + ^ 0, 

£/ien m defined via (|A.12[) . toztft. (1 — Co) = (1 — co)^ 1 and /i = (1 — co) _1 /2, solves 
the Riemann-Hilbert problem (jA.ll) and /1 = m±b± . 

Hence we have a formula for the solution of our Riemann-Hilbert problem m(z) 
in terms of (l—C w )~ 1 mo and this clearly raises the question of bounded invertibility 
of I — C w . This follows from Fredholm theory (cf. e.g. [38]): 

Lemma A. 6 ([22, 25 ). Assume Hvvothesis \A.1[ Then the operator I — C w is 
Fredholm of index zero, 

(A.15) fnd(I - C w ) = 0. 

By the Fredholm alternative, it follows that to show the bounded invertibility of 
I — C w we only need to show that ker(I — C w ) = 0. The latter being equivalent to 
unique solvability of the corresponding vanishing Riemann-Hilbert problem. 
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Corollary A. 7. Assume Hvpothesis \A.l[ 

A unique solution of the Riemann-Hilbert problem (|A.1|) exists if and only if the 
corresponding vanishing Riemann-Hilbert problem, where the normalization condi- 
tion is given by m(p 00 ) = (0 0), has at most one solution. 

We are interested in comparing two Riemann-Hilbert problems associated with 
respective jumps wq and w with \\w — Wo\\oo small, where 

(A.16) IHloo = || || + ||W-|| L oo (E) . 

For such a situation we have the following result: 

Theorem A. 8 f |24jb Assume that for some data Wq the operator 

(A.17) l-C w t : L 2 (£) -> L 2 (£) 

has a bounded inverse, where the bound is independent oft. 
Furthermore, assume w satisfies 

(A.18) - toJJIoo < 

for some function a(t) — > as t — > oo. Then (I — C w t)~ l : L 2 (S) — > L 2 (T,) also 
exists for sufficiently large t and the associated solutions of the Riemann-Hilbert 
problems (jA.ip only differ by 0{a(t)). 

Acknowledgments. We want to thank Ira Egorova for most helpful discussions 
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uscript. 
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